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We examine the impact of a non-minimal coupling of the inflaton to the Ricci scalar, |£7?<^ 2 , on 
the inflationary predictions. Such a non-minimal coupling is expected to be present in the infla¬ 
ton Lagrangian on fairly general grounds. As a case study, we focus on the simplest inflationary 
model governed by the potential V oc (j > 2 , using the latest combined 2015 analysis of Planck and 
BICEP2/AecA; Array. We find that the presence of a coupling £ is favoured at a significance of 
99% CL, assuming that nature has chosen the potential V oc (j> 2 to generate the primordial per¬ 
turbations and a number of e-foldings N = 60. Within the context of the same scenario, we find 
that the value of £ is different from zero at the 2<r level. When considering the cross-correlation 
polarization spectra from BICEP2/AecA; Array and Planck, a value of r = 0.038 /q'q 3 o * s predicted 
in this particular non-minimally coupled scenario. Future cosmological observations may therefore 
test these values of r and verify or falsify the non-minimally coupled model explored here. 

PACS numbers: 98.70.Vc, 98.80.Cq, 98.80.Bp 


I. INTRODUCTION 

Inflation provides the most theoretically attractive and 
observationally successful cosmological scenario able to 
generate the initial conditions of our universe, while solv¬ 
ing the standard cosmological problems. Despite this re¬ 
markable success, the inflationary paradigm is still lack¬ 
ing firm observational confirmation. The picture that 
emerges from the latest data from Planck, including also 
the joint analysis of R-mode polarization measurements 
from the BICEP2 collaboration [1-4], is compatible with 
the inflationary paradigm. According to these observa¬ 
tions, structure grows from Gaussian and adiabatic pri¬ 
mordial perturbations. From the theoretical viewpoint, 
this picture is usually understood as the dynamics of a 
single new scalar degree of freedom, the inflaton, min¬ 
imally coupled to Einstein gravity. However, the infla¬ 
ton (j> is expected to have a non-minimal coupling to the 
Ricci scalar through the operator f^Rcj) 2 , where £ is a 
dimensionless coupling. Indeed, successful reheating re¬ 
quires that the inflaton is coupled to the light degrees 
of freedom. Such couplings, though weak, will induce 
a non-trivial running for £. Thus, even starting from 
a vanishing value of £ (away from the conformal fixed 
point £ = —1/6) at some energy scale, a non-trivial non- 
minimal coupling will be generated radiatively at some 
other scale (see e.g. Ref. [5]). Therefore, it is important 
to study the impact of such a coupling on the inflationary 
predictions, especially in view of the latest Planck 2015 
data. 

Generically, for successful inflation, the inflaton should 
be very weakly coupled 1 * . It follows that the magnitude 


1 This requirement is also dictated by the non-detection of large 


of £ is expected to be small. Yet, even with such a sup¬ 
pressed coupling, the inflationary predictions are signifi¬ 
cantly altered [6-16]. For instance, and as we will see, a 
small and positive £ can enlarge considerably the space of 
phenomenologically acceptable scenarios (see also [17]). 
In this paper, we will focus on the simplest inflation¬ 
ary scenario with a potential V oc f> 2 [18], and a non¬ 
zero non-minimal coupling. According to the very recent 
Planck 2015 full mission results, the minimally-coupled 
version of this scenario (be. £ = 0) is ruled out at more 
than 99% confidence level [2, 4], for 50 e-folds of infla¬ 
tion. Nevertheless, the N = 60 case is only moderately 
disfavoured at 95% CL. Thus, before discarding it def¬ 
initely from the range of theoretical possibilities, it is 
worthwhile to explore this scenario in all generality (con¬ 
sidering as well different posibilities for the number of 
e-folds), given that, as explained earlier, the presence of 
non-minimal couplings in the inflaton Lagrangian is quite 
generic. 


II. NON-MINIMALLY COUPLED INFLATON 


The dynamics of a non-minimally coupled scalar field 
<j) with a potential U(</>) is governed, in the Jordan frame, 
by the following action 2 




R+ 


(1) 


primordial non-Gaussianities [3] and the soft breaking of the shift 
symmetry (f> — > <j> + c, necessary to protect the flatness of the 
potential. 

2 As usual, Mp = 1/y/SnGjy — 2.43 x 10 18 GeV is the reduced 

Planck mass. 
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where indices are contracted with the metric g defined 
as ds 2 = — df 2 + a 2 (t ) da; 2 . Inflation can be conveniently 
studied in the Einstein frame , after performing a confor¬ 
mal transformation g E v = fl(4>)g with f1 = l+££> 2 /M|, 
and canonically-normalizing the scalar field. Up to a to¬ 
tal derivative, the action takes the familiar form 

S = j d 4 a: y/~g E Re ~ ^e 9 ^ 9 ^ ~ v [ ( t ) (<?)] j . 

( 2 ) 

where now tp is the canonically-normalized inflaton, re¬ 
lated to the the original non-minimally coupled scalar 
field </> through 
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1 

n 



( 3 ) 


In terms of the original scalar field q i, the physical poten¬ 
tial takes the simple form 

U^(</>)] = U(^)/C! 2 (^). (4) 
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In the following, as previously stated, we shall focus 
on the simplest inflationary model. A generalization 
to other interesting inflationary scenarios, as for in¬ 
stance, the Higgs inflation model [19], will be carried 
out elsewhere [20]. The simplest scenario is given by 
the quadratic potential U(cf>) = \ m 2 q i> 2 , with a non¬ 
vanishing coupling £. In order to derive the primordial 
scalar and tensor perturbation spectra within the non- 
minimally coupled <f 2 theory, we shall make use of the 
slow-roll parameters' 5 : 






&R = 5) 


It is straightforward to derive the expressions for the 
spectral index of the primordial scalar perturbations 
n s = 1 + 2g — 6e, its running a = d n s / dlnfc = 
—24e 2 +16e?7—2 ^sr, and the tensor-to-scalar ratio r = 16e 
from the above slow-roll parameters 4 . 

Within the slow-roll approximation, one can easily 
solve numerically the inflationary dynamics governed by 
the action Eq. (2). The number of e-folds is given by 


1 

r* 1 
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3 Here, we use the notation £;s 11 (V 9 ) t° refer to the usual slow-roll 
parameter £, in order to avoid confusion with the non-minimal 
coupling to gravity £. 

4 Notice that the expressions for both n s and r are first-order in 

slow-roll, while a involves second order slow-roll terms. However, 

we have checked numerically that such second order corrections 

in slow-roll leave unchanged the constraints on the inflationary 

observables ( n s ,r ). Therefore, higher order slow-roll corrections 

can be safely neglected. 


FIG. 1. Theoretical predictions for the chaotic model V oc (j > 2 
with a non-minimal coupling £ in the (n s , r) plane for N = 50 
and N = 60. The red circles represent the £ = 0 case, corre¬ 
sponding to the usual predictions of the chaotic inflationary 
scenario. We show as well the 68% and 95% confidence level 
regions arising from the usual analyses in the (n s , r) plane 
using the various data combinations considered here. 


where e(<^>) = [V 7 (<j>)/V. 

The inflationary theoretical predictions for the N = 50 
and N = 60 cases are depicted in Fig. 1, in the (n s ,r) 
plane, for both positive and negative values of the cou¬ 
pling £. The case of £ = 0 corresponds to the usual 
predictions of the chaotic inflationary scenario, with 
n s = 1 - 2/N ~ 0.967 (n s ~ 0.96) and r = 8/N ~ 0.13 
(r ~ 0.16) for N = 60 (N = 50), and it is represented 
by red circles. Notice that negative values of £ lead to a 
larger tensor-to-scalar ratio. Positive values of £, on the 
other hand, will reduce the tensor contribution, while 
also pushing n s significantly below scale invariance as £ 
increases. For instance, for £ > 0.002 and N = 60, the 
scalar spectral index will always be smaller than the ob- 
servationally preferred value n s — 0.96. 

The predicted running of the spectral index a is shown 
in Fig. 2 as a function of the non-minimal coupling £. 
In general, negative (positive) values of £ lead to posi¬ 
tive (negative) values of the running. Although the large 
positive values of the running shown in Fig. 2 are com¬ 
patible with the recent Planck 2015 constraints [2], a = 
—0.0065 ± 0.0076, they are nevertheless associated with 
values of the tensor-to-scalar ratio r > 0.5, which are ex¬ 
cluded observationally. The red circle in Fig. 2 refers to 
the £ = 0 case, corresponding to a = —2/N 2 ~ —0.00056 
for N = 60. 















3 


III. OBSERVATIONAL CONSTRAINTS ON £ IN 

THE QUADRATIC INFLATIONARY MODEL 

In this paper, we restrict our numerical fits to Cosmic 
Microwave Background (CMB) measurements. The in¬ 
clusion of external data sets, such as Baryon Acoustic Os¬ 
cillation measurements, or a Hubble constant prior from 
the HST team will not affect the constraints presented 
in the following. Our data sets are the Planck tempera¬ 
ture data (hereafter TT) [21-23], together with the low-£ 
WMAP 9-year polarization likelihood, that includes mul¬ 
tipoles up to £ = 23, see Ref. [24] (hereafter WP), and 
the recent multi-component likelihood of the joint analy¬ 
sis of BICEP2/Aecfc Array and Planck polarization maps 
(hereafter BKP ), following the data selection and fore¬ 
ground parameters of the fiducial analysis presented in 
Ref. [1] However, variations of this fiducial model will 
not change significantly the results presented here. 

These data sets are combined to constrain the cosmo¬ 
logical model explored here, and described by the param¬ 
eters 5 6 * * * : 

{w b ,w c ,0 s ,T,log[lO lo A s ],£} (7) 

In Table I, we summarize the definition as well as the 
priors on these parameters. We use the Boltzmann code 
CAMB [25] and the cosmological parameters are extracted 
from the data described above by means of a Monte Carlo 
Markov Chain (MCMC) analysis based on the most re¬ 
cent version of cosmomc [26]. The constraints obtained 
on the non-minimal coupling £ are then translated into 
bounds on the usual inflationary parameters n B , r and a. 

Table II shows the 95% CL constraints on the param¬ 
eter £ as well as on the derived inflationary parameters 


Parameter 

Physical Meaning 

Prior 

ujb = Qb/i 2 

Present baryon density. 

0.005 -»• 0.1 

uj c = Q c h 2 

Present Cold dark matter density. 

0.01 0.99 

e s 

r B /D A (zdec) a - 

0.5 10 

T 

Reionization optical depth. 

0.01 -* 0.8 

ln(10 10 A fl ) 

Primordial scalar amplitude. 

2.7 -A 4 

£ 

Non-minimal coupling. 

-0.002 ->• 0.0065 


a The parameter © s is the ratio between the sound horizon r 3 
and the angular diameter distance Da (zd ec ) a t decoupling Zdec- 


TABLE I. Uniform priors for the cosmological parameters 
considered in the present analysis. 


5 This fiducial analysis assumes a tensor spectral index nj = 0, 
the BB bandpowers of BICEP2 /Keck Array and the 217 and 
353 GHz bands of Planck, in the multipole range 20 < t < 200. 

6 Notice that the inflationary cosmology under study contains less 

parameters than the standard ACDM picture, as once the non- 

minimal coupling £ is fixed, n s , r and a are fully determined, 

and are thus derived, parameters. 





FIG. 2. The running a as a function of the non-minimal 
coupling £. The red circle represents the minimal coupling 
case £ = 0. 


n s , r and the running a arising from our numerical analy¬ 
ses using the two CMB data combinations used here and 
assuming that n s and r are univocally determined by £ 
(for a fixed number of e-folds N , that we consider to be 
either 60 or 50). For N = 60, the preferred value of the 
non-minimal coupling £ from Planck TT plus WP mea¬ 
surements is positive and slightly larger than the mean 
value obtained when the cross-correlated polarized maps 
from BICEP2/AecA: and Planck (BKP) experiments are 
included in the numerical analyses. This preference for 
a slightly larger £ (and consequently, smaller r) is clear 
from the one-dimensional posterior probability distribu¬ 
tion of £ shown in the left panel of Fig. 4. The mean 
value of £ = 0.0028 obtained from Planck TT plus WP 
data is translated into a 95% CL constraint of the tensor- 
to-scalar-ratio r = 0.038+ 0 ; 031 , as can be seen from the 
right panel of Fig. 4. When considering BICEP2/ Keck 
and Planck cross-spectra polarization data, the former 
constraint on the tensor-to-scalar ratio is very similar to 
the one quoted above. Concerning the running of the 
spectral index, the two data combinations seem to have 
a preference for a small negative running a = —0.0005, 
associated to small values of [£[, as shown in Fig. 2. 

Let us now comment on the sensitivity of our con¬ 
straints to changes in the number of e-folds N. Setting 
N = 50 leads to different, though almost insignificant, 
changes in the constraints obtained using the two CMB 
data sets. The theoretically allowed regions in the (n s ,r) 
plane as a function of £ for N = 50 are indeed slightly 
different from those corresponding to the N = 60 case, 
see Fig. 1. The net result is a smaller (larger) values 
of n s (r) than in the N = 60 case. The BICEP2/Aecfc 
and Planck cross-spectra polarization data yield a value 
r = 0.053 ^o'q 37 for the tensor-to-scalar ratio in the con¬ 
text of the non-minimally coupled </> 2 model. On the 
other hand, the resulting central value for the scalar spec¬ 
tra index is only half a a away (towards smaller values) 
from the corresponding one for N = 60, as expected from 
the theoretical predictions illustrated in Fig. 1. 












A< p/Mp 
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TABLE II. Inflationary constraints in the context of non-minimally coupled chaotic potential <j> 2 -. The upper block 
of the table refers to the 95% CL limits on the non-minimal coupling £ (the parameter varied in the MCMC analyses) from the 
two possible CMB data combinations used in this study, for both N = 60 and N = 50. The lower block of the table contains 
the 95% CL derived ranges of the inflationary parameters n s , r and a from the limits of £ illustrated above, in the context of 
the non-minimally coupled chaotic potential <f> 2 , for both N = 60 and N = 50. 


Figure 1 shows the 68% and 95% CL allowed regions 
in the ( n s ,r ) plane resulting from our MCMC analyses 
to Planck TT plus WP data and to the combined BKP 
in the usual ( n s ,r ) plane, together with the theoretical 
predictions for N = 50 and N = 60 for the non-minimally 
coupled <f) 2 scenario. 

To address the question of whether or not a non- 
minimal coupling £ is favoured by current CMB data, we 
compare the x 2 test statistics function for the </> 2 model 
in its minimally and non-minimally coupled versions for 
N = 60, albeit very similar results are obtained for 



• V (X <p, i = 0, N = 60 




11 12 13 14 15 16 

A (j)/Mp 



0.08 0.16 0.24 0.32 0.40 0.48 0.50 0.64 


r 

FIG. 3. Excursion of the canonically-normalized inflaton 
versus the one of the original scalar <f>. The magnitude of 
the tensor-to-scalar ratio is encoded in the curve through the 
color bar. Notice that, in both frames, large r correlates with 
large excursions. 


N = 50. The x 2 for the case of Planck TT plus WP data, 
evaluated at the best-fit-point of the </> 2 model minimally 
coupled to gravity is x 2 [£ = 0] = 9812.8. On the other 
hand, the non-minimally coupled version has a lower x 2 
value at the best-fit-point due to the extra parameter £ 
introduced in the model, with x 2 [£ ^ 0] = 9806.8. The 
difference between these two x 2 values is Ax 2 = 6, which, 
for a distribution of one degree of freedom, has a p -value 
of 0.014, and is considered as statistically significant. For 
the case of the combined BKP likelihood, the difference 
between the test-statistics for the minimally coupled and 
non-minimally coupled <j) 2 models is Ax 2 = 10, which, 
for one degree of freedom, has a p -value of 0.0016, and is 
considered as very statistically significant. Therefore, ac¬ 
cording to the most recent CMB data, the presence of a 
non-minimal coupling £ within the (j) 2 model is favoured 
at a significance equal or larger than ~ 99% CL. 

Let us now turn to future constraints on £. Future 
observations, as those expected from PIXIE [27], Eu¬ 
clid [28], COrE [29] and PRISM [30], could be able to 
reach an accuracy of ay = cr„ s _i = 10 -3 . With such 
precision, one could hope to test deviations from the 
quadratic potential [31], as the one studied here, by con¬ 
structing quantities independent of N, up to subleading 
0(1/N 3 ) corrections. It is straighforward to get for our 
case, 

n s ~ 1 + | =~ 20£, (8) 

at leading order both in slow-roll and £. If it turns out 
that nature had chosen a very small value of r, future 
constraints on £ would be as strong as £ < ICC 4 ; one 
order of magnitude stronger than the ones obtained in 
this analysis. Concerning the running a, it is interesting 
to note that futuristic observations like SPHEREx [32] 
with a forcasted error of cr a = 10~ 3 , will be able to falsify 
the present scenario. 

Finally, it is also interesting to explore the impact of 
the non-minimal coupling on the inflaton excursion. It is 
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FIG. 4. The left (right) panel shows the one-dimensional posterior probability distributions of the non-minimal coupling £ (the 
tensor-to-scalar ratio r) in the context of a non-minimally coupled chaotic potential rj> 2 , with r a function of £, and therefore, 
a prediction within the model. 


well-know that large values of the tensor-to-scalar ratio 
r, as those found by previous BICEP2 measurements ' 
[33, 34] yield large inflaton excursions </ M P [35-38], 
which are hard to understand in the context of a con¬ 
sistent effective held theory. In particular, successful in¬ 
flation requires that higher order non-renormalizable op¬ 
erators, which are expected to be naturally present in 
the inhationary potential, are sufficiently suppressed. A 
number of phenomenological studies have recently been 
devoted to address this problem [39-42]. In Fig. 3, we 
plot the excursion of both </> and p, together with the 
corresponding tensor-to-scalar ratio r. It is clear that 
the excursion of the canonically-normalized inhaton p is 
lowered for positive values of £ i.e. Ap < A</>. How¬ 
ever, this decrease is rather mild and the excursion still 
takes on super-Planckian values for the phenomenologi¬ 
cally acceptable values of £. Conversely, negative values 
of £ lead to an increase of the excursion of p. Figure 3 
also shows that super-Planckian values of both <f> and 
p are still associated with large values of the tensor-to- 
scalar ratio r, in agreement with the Lyth bound [35]. 
Thus, once a small non-zero and positive value of the 
coupling £ is turned on, both the inflaton excursion and 
r are slightly lowered, but without alleviating completely 
the super-Planckian excursion problem. 

IV. CONCLUSIONS 

A small, non-minimal coupling ^R<f> 2 is expected to 
be present in the inflaton Lagrangian, and modifies the 


' The joint BKP analysis finds however no evidence for primordial 


inflationary predictions in an interesting way. Focusing 
on the simplest quadratic potential scenario, and using 
the very recent joint analysis of BICEP2/ATecfc Array and 
Planck polarization maps, we found that a small, posi¬ 
tive value of the coupling £ is favoured at the 2er level, 
assuming that nature has chosen the <fi 2 scenario for the 
generation of primordial perturbations. If only Planck 
TT plus WP data are used in the analyses, the signif¬ 
icance is milder. These conclusions have been obtained 
for a number of e-foldings within the N = 50 — 60 range. 
It would be interesting to see if upcoming U-modes mea¬ 
surements can reinforce or weaken the statistical signifi¬ 
cance of these findings. In particular, it would be crucial 
to discriminate between the presence of a non-minimal 
coupling in the theory and other departures from the 
quadratic approximation. 
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6 


[1] P. A. R. Ade et al. [BICEP2 and Planck Collaborations], 
“A Joint Analysis of BICEP2/Keck Array and Planck 
Data,” [arXiv:1502.00612 [astro-ph.CO]]. 

[2] [Planck Collaboration], “Planck 2015 results. XIII. Cos¬ 
mological parameters,” arXiv: 1502.01589 [astro-ph.CO]. 

[3] P. A. R. Ade et al. [Planck Collaboration], Planck 
2015 results. XVII. “Constraints on primordial non- 
Gaussianity,” arXiv:1502.01592 [astro-ph.CO]. 

[4] P. A. R. Ade et al. [Planck Collaboration], “Planck 2015. 
XX. Constraints on inflation,” arXiv: 1502.02114 [astro- 
ph.CO]. 

[5] I. L. Buchbinder, S. D. Odintsov and I. L. Shapiro, Effec¬ 
tive action in quantum gravity, Bristol, UK: IOP (1992) 
413 pp. 

[6] D. S. Salopek, J. R. Bond and J. M. Bardeen, Phys. Rev. 
D 40, 1753 (1989). 

[7] T. Futamase and K. i. Maeda, Phys. Rev. D 39 (1989) 
399. 

[8] R. Fakir and W. G. Unruh, Phys. Rev. D 41, 1783 (1990). 

[9] D. I. Kaiser, Phys. Rev. D 52, 4295 (1995) [astro- 
ph/9408044], 

[10] E. Komatsu and T. Futamase, Phys. Rev. D 59, 064029 
(1999) [astro-ph/9901127]. 

[11] M. P. Hertzberg, JHEP 1011, 023 (2010) 

[arXiv: 1002.2995 [hep-ph]]. 

[12] N. Okada, M. U. Rehman and Q. Shafi, Phys. Rev. D 82 
(2010) 043502 [arXiv:1005.5161 [hep-ph]]. 

[13] A. Linde, M. Noorbala and A. Westphal, JCAP 1103, 
013 (2011) [arXiv:1101.2652 [hep-th]]. 

[14] D. I. Kaiser and E. I. Sfakianakis, Phys. Rev. Lett. 112 
(2014) 1, 011302 [arXiv: 1304.0363 [astro-ph.CO]]. 

[15] T. Chiba and K. Kohri, PTEP 2015, no. 2, 023E01 
[arXiv: 1411.7104 [astro-ph.CO]]. 

[16] C. Pallis and Q. Shafi, JCAP 1503, no. 03, 023 (2015) 
[arXiv: 1412.3757 [hep-ph]]. 

[17] S. Tsujikawa, J. Ohashi, S. Kuroyanagi and A. De Fe¬ 
lice, Phys. Rev. D 88 (2013) 2, 023529 [arXiv: 1305.3044 
[astro-ph.CO]]. 

[18] A. D. Linde, Phys. Lett. B 129 (1983) 177. 

[19] F. L. Bezrukov and M. Shaposhnikov, Phys. Lett. B 659, 
703 (2008) [arXiv:0710.3755 [hep-th]]. 

[20] L. Boubekeur, E. Giusarma, O. Mena and H. Ramirez, 
In preparation. 

[21] P. A. R. Ade et al. [Planck Collaboration], 

arXiv: 1303.5062 [astro-ph.CO]. 

[22] P. A. R. Ade et al. [Planck Collaboration], 

arXiv:1303.5075 [astro-ph.CO]. 


[23] P. A. R. Ade et al. [Planck Collaboration], 
arXiv: 1303.5077 [astro-ph.CO]. 

[24] C. L. Bennett et al. [WMAP Collaboration], Astrophys. 
J. Suppl. 208, 20 (2013) [arXiv: 1212.5225 [astro-ph.CO]]. 

[25] A. Lewis, A. Challinor and A. Lasenby, Astrophys. J. 
538, 473 (2000) [arXiv:astro-ph/9911177]. 

[26] A. Lewis and S. Bridle, Phys. Rev. D 66, 103511 (2002) 
[arXiv:astro-ph/0205436]. 

[27] A. Kogut, D. J. Fixsen, D. T. Chuss, J. Dotson, E. Dwek, 
M. Halpern, G. F. Hinshaw and S. M. Meyer et al, JCAP 
1107 (2011) 025 [arXiv: 1105.2044 [astro-ph.CO]]. 

[28] R. Laureijs et al. [EUCLID Collaboration], “Euclid Def¬ 
inition Study Report,” arXiv:1110.3193 [astro-ph.CO]. 

[29] F. R. Bouchet et al. [COrE Collaboration], “COrE (Cos¬ 
mic Origins Explorer) A White Paper,” arXiv: 1102.2181 
[astro-ph.CO]. 

[30] P. Andre et al. [PRISM Collaboration], “PRISM (Po¬ 
larized Radiation Imaging and Spectroscopy Mission): 
A White Paper on the Ultimate Polarimetric Spectro- 
Imaging of the Microwave and Far-Infrared Sky,” 
arXiv: 1306.2259 [astro-ph.CO]. 

[31] P. Creminelli, D. Lopez Nacir, M. Simonovic, G. Tre- 
visan and M. Zaldarriaga, Phys. Rev. Lett. 112 (2014) 
24, 241303 [arXiv: 1404.1065 [astro-ph.CO]]. 

[32] O. Dore, J. Bock, P. Capak, R. de Putter, T. Eifler, 
C. Hirata, P. Korngut and E. Krause et al., “SPHEREx: 
An All-Sky Spectral Survey,” arXiv: 1412.4872 [astro- 
ph.CO]. 

[33] P. A. R. Ade et al. [BICEP2 Collaboration], Phys. 
Rev. Lett. 112 (2014) 241101. [arXiv:1403.3985 [astro- 
ph.CO]]. 

[34] P. A. R. Ade et al. [BICEP2 Collaboration], Astrophys. 
J. 792 (2014) 62 [arXiv: 1403.4302 [astro-ph.CO]]. 

[35] D. H. Lyth, Phys. Rev. Lett. 78 (1997) 1861 [hep- 
ph/9606387]. 

[36] L. Boubekeur, Phys. Rev. D 87 (2013) 6, 061301 
[arXiv: 1208.0210 [astro-ph.CO]]. 

[37] G. Efstathiou and K. J. Mack, JCAP 0505 (2005) 008 
[astro-ph/0503360]. 

[38] L. Verde, H. Peiris and R. Jimenez, JCAP 0601 (2006) 
019 [astro-ph/0506036]. 

[39] J. Garcia-Bellido, D. Roest, M. Scalisi and I. Zavala, 
JCAP 1409, 006 (2014) [arXiv:1405.7399 [hep-th]]. 

[40] G. Barenboim and O. Vives, arXiv: 1405.6498 [hep-ph]. 

[41] J. Garcia-Bellido, D. Roest, M. Scalisi and I. Zavala, 
Phys. Rev. D 90, no. 12, 123539 (2014) [arXiv: 1408.6839 
[hep-th]]. 

[42] L. Boubekeur, E. Giusarma, O. Mena and H. Ranhrez, 
arXiv: 1411.7237 [astro-ph.CO]. 



